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Abstract 

O 

We propose to use rich informations on the pp,pp total cross sections o~tot below 
iV(~ 10)GeV in order to predict the total cross section and p ratio at very high 
energies. Using the FESR as a constraint for high energy parameters, we search for 
the simultaneous best fit to the data points of o" t ot and p ratio up to some energy 
(e.g., ISR, Tevatron) to determine the high-energy parameters. We then predict 
ctot and p in the LHC and high-energy cosmic-ray regions. Using the data up to 
i/s = 1.8TeV(Tevatron), we predict a^ t and p pp at the LHC energy(y / s = 14TeV) 
as 106.3 ± 5.1 syst ± 2.4 stat mb and 0.126 ± 0.007 syst ± 0.004 stat , respectively. The 
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predicted values of o~tot m terms of the same parameters are in good agreement 
with the cosmic-ray experimental data sample up to Pi a b ~ 10 8 ~ 9 GeV by Block, 
Halzen, and Stanov. 
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Recently [1], we have proposed to use rich informations on rep total cross sec- 
tions below N(~10 GeV) in addition to high-energy data to discriminate 
whether these cross sections increase like log v or log 2 v at high energies [2]. 
The FESR which was derived in the spirit of the P' sum rule [3] as well as the 
n — 1 moment FESR([4], [5]) have been required to constrain the high-energy 
parameters. We then searched for the best fit of <j[^ above 70GeV in terms 
of high-energy parameters constrained by the two FESR. We then arrived at 
the conclusion that our analysis prefers the log 2 !/ behaviours consistent with 
the Froissart-Martin unitarity bound [6]. 

As for the pp and pp total cross sections, there are a lot of data including 
cosmic-ray data up to yfs ~ several times of 10 4 GeV compared with data up 
to yfs ~30GeV for ttN scattering. Therefore, it is very valuable if one could 
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investigate the high-energy behaviours at LHC and cosmic- ray regions [8] using 
the similar approach as ref. [1]. 

(The purpose of this Letter): The purpose of this Letter is to predict a[„t ■, 
the pp, pp total cross sections and p^ + \ the ratio of the real to imaginary 
part of the forward scattering amplitude at the LHC and the higher-energy 
cosmic-ray regions, using the experimental data for <j[^ and p^ for 70GeV< 
Piab < Piarge as inputs. We first choose Piarge = 2100GeV corresponding to ISR 
region (-^/s ~ 60GeV). Secondly we choose Pi arg e = 2 x 10 6 GeV corresponding 
to the Tevatron collider (a/s ~ 2TeV). In a recent paper, Block and Halzen[7] 
emphasized the importance of p for the evidence for saturation of the Froissart- 
Martin bound[6]. We also use the p ratio as input data in addition to FESR 
as a constraint. We searched for the simultaneous best fit of and p^ in 
terms of high-energy parameters co,Ci,C2 and j3p> constrained by the FESR. 
It turns out that the prediction of a[„t agrees with pp experimental data at 
these cosmic-ray energy regions[8,22] within errors in the first case ( ISR ). It 
has to be noted that the energy range of predicted <7tot\ p^ is several orders 
of magnitude larger than the energy region of 0tot\ input (see Fig. 1). If 
we use data up to Tevatron (the second case), the situation is much improved, 
although there are some systematic uncertainties coming from the data at 
y/s = 1.8TeV (see Fig. 2). 

FESR(l): Firstly we derive the FESR in the spirit of the P' sum rule [3]. Let 
us consider the crossing-even forward scattering amplitude defined by 

FM[V) = f*» + fW) with /m F(+ , M = *fgM . (1) 

We also assume 

Im F i+) (u) = Im R(u) + Im F P ,(u) 

V ( , V , o V \ PP' ( V \ ap ' 

at high energies (v > N). We have defined the functions R{v) and Fpi{v) by 
replacing p by M in Eq. (3) of ref.[l]. Here, M is the proton( anti-proton) 
mass and u, k are the incident proton(anti-proton) energy, momentum in the 
laboratory system, respectively. 

Since the amplitude is crossing-even, we have 
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+ H log2£ AT + log2 s)} , 

Pp> f(e- in v/M) a p' + (u/M) a P r 

Fpr(v) = — — : 

v ' M \ sni7rQ!jv 
and subsequently obtain 
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Re f p ,(v) = -^- — . 
v ' M \MJ 
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substituting api — \ in Eq. (4). Let us define 

F i+ \v) = F i+ \v) - R(v) - F P ,(v) ~ (a(0) < 0) 
Using the similar technique to ref.[l], we obtain 
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where TV = ViV 2 - M 2 ~ JV. Let us call Eq. (8) as the FESR(l). 
(FESR(2)): The second FESR corresponding to n = 1 [5] is: 
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vim R[y)dv + ^ vim Fp>(v)dv 



We call Eq. (9) as the FESR(2) which we use in our analysis. 

(The ratio): The ratio, the ratio of the real to imaginary part 
F^ + \v) is obtained from Eqs. (2), (5) and (6) as 
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Re F(+) {y) Re R{u) + Re F P , (u) 
Im F(+)(z/) ~~ 7m + Im F P >{u) 



(ci + 2c 2 log^) - ^ (£) 



(10) 




(General approach): The FESR(l)(Eq. (8)) has some problem, i.e., there 
are the so-called unphysical regions coming from boson poles below the pp 
threshold. So, the contributions from unphysical regions of the first term of the 
right-hand side of Eq. (8) have to be calculated. Reliable estimates, however, 
are difficult. Therefore, we will not adopt the FESR(l). 

On the other hand, contributions from the unphysical regions to the first term 
of the left-hand side of FESR(2)(Eq. (9)) can be estimated to be an order of 
0.1% compared with the second term. 1 Thus, it can easily be neglected. 

Therefore, the FESR(2)(Eq. (9)), the formula of ^(Eqs. (1) and (2)) and 
the p( + "> ratio (Eq. (10)) are our starting points. Armed with the FESR(2), 
we express high-energy parameters Co, ci, c 2 , f3p> in terms of the integral of 
total cross sections up to N. Using this FESR(2) as a constraint for f3p> = 
/3p/(co, Ci, c 2 ), the number of independent parameters is three. We then search 
for the simultaneous best fit to the data points of cr^(k) and p^ + \k) for 
70GeV< k < Piarge to determine the values of c ,Ci,c 2 giving the least \ 2 ■ 
We thus predict the <j tot and in LHC energy and high-energy cosmic-ray 
regions. 

( Data ): We use rich data[9] of a pp and a pp to evaluate the relevant integrals 
of cross sections appearing in FESR(2). We connect the each data point 2 of 

1 The average of the imaginary part from boson resonances below the pp thresh- 
old is the smooth extrapolation of the t-channel qqqq exchange contributions 
from high energy to v < M due to FESR duality [4, 5]. Since Im F^~{v) < 
Im F(+V), ^ vim F^g(v)dv < ^ vim F^\v)dv = %Im f pp (v)dv ~ 

Im f pp \ k=0 ^ 3.2GeV < ^ f" k 2 a[^{k)dk = 3403±20GeV, where we use the 
experimental value, ^<J PP t ~14.4GeV _1 in k < 0.3GeV. So, resonance contributions 
to the first term of Eq. (9) is less than 0.1% of the second term. 
Besides boson resonances, there may be additional contributions from multi-pion 
contributions below pp threshold. In the pp annihilation, pp — > -kit could give com- 
parable contributions with p-meson, but multi-pion contributions are suppressed 
due to the phase volume effects. Therefore, the first term of Eq. (9) will still be 
negligible even if the above contributions are included. 



2 We take the error Ay for each data point y as Ay = y (Ay) 2 tat + (Ay) 2 yst . 

When several data points, denoted yi with error Ay« (i = l,---,n), are listed 
at the same value of k, these points are replaced by y with Ay, given by y = 
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k 2 <Jtot an d ^ 2 c"tot with the next point by a straight line in order, from k — to 
k = N, and regard the area of this polygonal line graph as the relevant integral 
in the region < k < N. The integral of k 2 a^t(k) is given by averaging those 
of k 2 aZ{k) and k 2 aZ(k). We have obtained 



f k 2 a ( t ^{k)dk = 3403 ± 20 GeV. (11) 
o 

for N = lOGeV (which corresponds to \/s = E cm = 4.54GeV). 3 The error of 
the integral, which is from the error of each data point, is very small (less than 
1%), and thus, we regard the central value as an exact one in the following 
analysis. 

When ofo t and of^ t data points are listed at the same value of k, we make the 
a iot(k) data point by averaging these values. Totally, 37 points are obtained 
in the energy region, 0.54GeV < k < 2100GeV. The data point of maximum 
value k = 2094.03GeV(v/s = 62.7GeV) comes from ISR[10]. There are 12 
points in the 70GeV< k < 2100GeV (11.5GeV < < 62.7GeV). There are 
no data reported in the wide range of 2100GeV < k < 1 x 10 5 GeV. There 
are 6 points[ll,12,15,14,13] of a^t in the Tevatron-collider energy region, 1 x 
10 5 GeV < k < 2 x 10 6 GeV. 

It is necessary to pay special attention to treat the data with the maxi- 
mum k = 1.7266 x 10 6 GeV(y / i = 1.8TeV) in this energy range, which comes 
from the three experiments E710[13]/E811[14] and CDF[15]. The former two 
experiments are mutually consistent and their averaged pp cross section is 
(jfo t = 72.0 ± 1.7mb, which deviates from the result of CDF experiment 
aZ = 80.03 ± 2.24mb. 

Again there are no data reported in the range 2 x 10 6 GeV < k < 2 x 10 7 GeV. 
There are 7 points of ofot with somewhat large errors, reported in the cosmic- 
ray energy region, 2 x 10 7 GeV <K5x 10 8 GeV (6TeV < ^ < 30TeV), 
coming from cosmic-ray experiments [8, 22]. Totally we obtain 25 points of 
in k > 70GeV. We have not included the cosmic-ray data in our analysis. 
Thus, 18 points of a^t are used in the analyses. 



[E i y i /(Ay i ) 2 ]/[E i V(A^) 2 ] and Ay = ^VE, V(Ay,) 2 ]- Then, the data points 
with Ay less than 3 mb are picked up. As a result, we obtain the 255(124) points for 
k 2 aZ t (k 2 aZ), giving the integrals (5.070 ± 0.034) x 10 4 ( (3.482 ± 0.037) x 10 4 )GeV 
in the region < k < N with N = lOGeV. 

3 The laboratory momentum are related to the CM energy squared s by s = 
2M(M + ^M 2 + P 2 ah ) and equivalently Pi ab = ^ yT - 4M 2 /s. Thus, at high 
energies E cm = y/s ~ xJ2MP Xoh . 
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The data of p pp (k)(= Re f pp (k)/Im f pp {k)) and p pp (k)(= Re f pp (k)/Im f pp (k)) 
are reported in ref.[9]. When both data points are listed at the same value 
of k, we can make the p {+) {k)(= Re F {+ \k)/Im F^ + \k)) data point. 4 We 
obtain 9 points of p( + ) in the energy region, 70GeV < k < 2100GeV. 5 
No data are reported in the range 2100GeV < k < 1 x 10 5 GeV. The two 
points of fP p are reported in Tevatron-collider energy region, 1 x 10 5 GeV < 
k < 2 x 10 6 GeV ( at k = 1.5597 x 10 5 GeV(v^ = 541GeV)[17] and k = 
1.7266 x 10 6 GeV( v / i =1.8TeV)[13] ). We regard these two points as the 
data. As a result, we obtain 11 points of pW up to Tevatron-collider energy 
region, 70GeV < k < 2 x 10 6 GeV. 

In the actual analyses, we use Re F(+) instead of p (+) (= Re F^/Im fW). 
The data points of Re (k) are made by multiplying p^ (k) by Im F^ (k) = 
^(°tot(k) + Otot(k))- The values of a{ p t and a PP t at the relevant values of k are 
obtained as follows: For k < 1500GeV, they are determined by the formula 
given in ref.[9](see the footnote 4). Two experimental values[12,13] of a pp in 
the Tevatron region are used. 

( Analysis ): As was explained in the general approach, both a[„t an d Re F^ 
data in 70GeV < k < Pi arg e are fitted simultaneously through the formula 
Eq. (2) and Eq. (10) with the FESR(2)(Eq. (9)) as a constraint. FESR(2) 
with Eq. (11) gives us 

8.87 = c + 2.04 Cl + 4.26c 2 + 0.367/3 P , , (12) 

which is used as a constraint of (3 pi = /3p'(c , c±, c 2 ), and the fitting is done by 
three parameters c ,Ci and c 2 . 

We have done for the following three cases: 

fit 1): The fit to the data up to ISR energy region, 70GeV < k < 2100GeV, 
which includes 12 points of a^t an d 7 points of p^ + \ 

fit 2): The fit to the data up to Tevatron-collider energy region, 70GeV< 
k < 2 x 10 6 GeV. For k = 1.7266 x 10 6 GeV(v^ = 1.8TeV), the E710/E811 
datum is used. There are 18 points of <j[^ and 9 points of p^ + \ 
fit 3): The same as fit 2, except for the CDF value at y/s = 1.8TeV used. 

( Results of the fit ): The results are shown in Fig. l(Fig. 2) for the fit l(fit 

4 Here the values of Im f pp (k) and Im f pp (k) at the relevant values of k are 
determined through the formula given in ref.[9], aPP/PP = Z + B\og 2 (s/s ) + 
Yi^i/sf' 1 ± Y 2 (s 1 /s)n with Yi,y 2 )=(35.45, 0.308, 42.53, 33.34)mb, 
( v ^o.V^)=( 5 - 38 ' 1 ) GeV and (m,%)=(0.458, 0.545). 

5 Here only the data point of maximum k = 1479GeV(\/i = 52.7GeV) is ob- 
tained by combining the ffP at k = U73A6GeV(^s = 52.6GeV) and at 
k = 1484.69GeV(Vs = 52.8GeV), reported in ref.[16]. The other 8 points are ob- 
tained by combining fP p and fP p with the same values of k. 
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Table 1 

The values of \ 2 f°r the fit f (fit up to ISR energy) and the fit 2 and fit 3 (fits 
up to Tevatron-collider energy). Np and N a (N p ) are the degree of freedom and the 
number of cr^(p^) data points in the fitted energy region. 





X 2 /N F 




X 2 p /N p 


fit 1 


10.6/15 


3.6/12 


7.0/7 


fit 2 


16.5/23 


8.1/18 


8.4/9 


fit 3 


15.9/23 


9.0/18 


6.9/9 



Table 2 



The best-fit values of parameters in the fit 1, fit 2 and fit 3. 





C2 


Cl 


CO 


0p> 


fit 1 


0.0411 ±0.0199 


-0.074 T 0.287 


5.92 ± 1.07 


7.96 =F 1.55 


fit 2 


0.0412 ±0.0041 


-0.076 T 0.069 


5.93 ±0.28 


7.95 T 0.44 


fit 3 


0.0484 ± 0.0043 


-0.181 T 0.071 


6.33 ±0.29 


7.37 T 0.45 



2 and fit 3). The x 2 /d-o.f are given in Table 1. The reduced x 2 an d the 
respective x 2 -values devided by the number of data points for and 
are less than or equal to unity. The fits are successful in all cases. There 
are some systematic differences between fit 2 and fit 3, which come from the 
experimental uncertainty of the data at y/s = 1.8TeV mentioned above. 

The best-fit values of the parameters are given in Table 2. Here the errors of 
one standard deviation are also given. 6 

(Predictions for and p^ at LHC and cosmic-ray energy region): By us- 
ing the values of parameters in Table 2, we can predict the and p( + ) in 
higher energy region, as are shown, respectively in (c) and (d) of Fig. 1 and 
2. The thin dot-dashed lines represent the one standard deviation. 

As is seen in (c) and (d) of Fig. 1, the fit 1 leads to the prediction of a^t 
and p( + ) with somewhat large errors in the Tevatron-collider energy region, 
although the best-fit curves are consistent with the present experimental data 
in this region. Furthermore, the predicted values of agree with pp exper- 
imental data at the cosmic-ray energy regions[8,22] within errors (see (a),(c) 
of Fig. 1). The best-fit curve gives x 2 /(number of data) to be 13.0/16, and 



6 The C2log 2 (^/M)-term in Eq. (2) is most relevant for predicting a[ ot in high 
energy region, and thus, the error estimation is done as follows: The is fixed 
with a value deviated a little from the best-fit value, and then the x 2_n t is done 
by two parameters c$ and c±. When the resulting \ 2 is larger than the least \ 2 
of the three-parameter fit by one, the corresponding values of parameters give one 
standard deviation. 
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(c) a[*t '■ High energy region (d) p^ + > 

Fig. 1. Predictions for o"W an d p(+) [ n terms of the fit 1. The fit is done for the 
data up to the ISR energy, in the region 70GeV< k < 2100GeV (11.5GeV < y/s < 
62.7GeV) which is shown by the arrow in each figure. Total cross section 
in (a) all energy region, versus logio-FLft/GeV, (b) low energy region (up to 
ISR energy), versus Pi ab /GeV and (c) high energy (Tevatron-collider, LHC 
and cosmic-ray energy) region, versus center of mass energy E cm in TeV unit, 
(d) gives the p^ + \— Re F^/Im F^) in high energy region, versus E cm in 
terms of TeV. The thin dot-dashed lines represent the one standard deviation. 

the prediction is successful. As was mentioned in the purpose of this Letter, 
it has to be noted that the energy range of predicted a[^t is several orders of 
magnitude larger than the energy region of the , input. If we use data 
up to Tevatron-collider energy region as in the fit 2 and fit 3, the situation is 
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mb 




(c) cr^t : High energy region (d) p^ + > 

Fig. 2. Predictions for and in terms of the fit 2(shown by green lines) 
and fit 3(shown by blue lines). The fit is done for the data up to Tevatron-collider 
energy, in the region 70GeV< k < 2 x 10 6 GeV(11.5GeV < y/s < 1.8TeV) which 
is shown by the arrow. For k = 1.7266 x 10 6 GeV(y / s = E cm = 1.8TeV), the aver- 
aged datum of E710[13]/E811[14], ofg. = 72.0 ± 1.7mb, is used in fit 2, while the 
a tot = 80.03 db 2.24mb of CDF[15] is used in fit 3. For each figure, see the caption 
in Fig.l. 

much improved (see (a),(c) of Fig. 2), although there is systematic uncertainty 
depending on the treatment of the data at \fs = 1.8TeV. 

The best-fit curve gives x 2 / (number of data) from cosmic-ray data, 1.3/7(1.0/7) 
for fit 2(fit 3). 
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Table 3 



The predictions of a^J and at LHC energy y/s = E cm = 

14TeV(P iafc =1.04xl0 8 GeV), and at a very high energy P lab = 5 • 10 20 eV 
( v / s=£' cm =967TeV.) in cosmic-ray region. 







(v^=14TeV) 




/O (+) (Piaf,=5-10 20 eV) 


fit 1 


103.8 ± 14.3mb 


0.122 ±£8SS 


188 ± 43mb 


0.099 


fit 2 


103.8 ± 2.3mb 


0.122 ±0.004 


189 ± 8mb 


0.100 ±0.003 


fit 3 


108.9 ± 2.4mb 


0.129 ±0.004 


204 ± 8mb 


0.104 ±0.003 



We can predict the values of a[tt an d a t LHC energy, v /i=£ , cm =14TeV 
and at very high energy of cosmic-ray region. The relevant energies are very 
high, and the and p^ can be regarded to be equal to the of<f t and p pp . 
The results are shown in Table 3. 

The prediction by the fit 1 in which data up to the ISR energy are used as input 
has somewhat large (fairly large) errors at LHC energy (at high energy of cosmic 
ray). By including the data up to the Tevatron collider, the prediction of fit 
2(using E710/E811 datum) is smaller than that of fit 3(using CDF datum). 
We regard the difference between the results of fit 2 and fit 3 as the systematic 
uncertainties of our predictions. As a result, we predict 

aZ = 106.3 ± 5.1 syst ± 2.4 stat mb, ff p = 0.126 ± 0.007 syst ± 0.004 stat (13) 

at LHC energy(y / i = E cm = 14TeV). We obtain fairly large systematic errors 
coming from the experimental unceratinty at ^/s = 1.8 TeV. 

The predicted central value of a^t is in good agreement with Block and 
Halzen[18] of? t = 107.4±1.2 mb, p pp = 0.132±0.001. In contrary to our results( 
see Fig. 2(a), (c)), however, their values are not affected so much about CDF, 
E710/E811 discrepancy. Our prediction has also to be compared with Cudell 

et al.[23] ergf t = 111.5 ± 1.2 syst -2.1^ mb, p pp = 0.1361 ± 0.0015 syst ±aoo25 stat , 
who's fitting techniques favours the CDF point at ^fs = 1.8 TeV. 

Finally we emphasize that precise measurements of both cr^t and pP v in the 
coming LHC experiments will resolve the FNAL discrepancy of ofo t ( Fig. 2(a), 
(c)). The LHC measurements would also clarify which is the best solution 
among the three high-energy cosmic-ray samples[20, 21,22]. 
Note added in proof: After completion of the hep-ph/0505058, we were in- 
formed that M. M. Block and F. Halzen[18] have also done the similar work 
based on the same spirit of duality using different method independently. We 
were also informed by M. J. Menon[19] about other cosmic-ray analyses by 
Gaisser et al.[20] and N. N. Nikolaev[21] besides M. M. Block et al.[22] which 
are used as cosmic-ray data in this Letter. 
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